Damping of electromagnetic waves due to electron-positron pair production 



I. INTRODUCTION 



S. S. Bulanov/'E A. M. Fedotov, 2 'Q and F. Pegoraro 3 '0 

institute of Theoretical and Experimental Physics, 117218 Moscow, Russia 
2 Moscow State Engineering Physics Institute, 115409 Moscow, Russia 
'Department of Physics, University of Pisa and INFM, Pisa, Italy 

The problem of the backreaction during the process of electron-positron pair production by a 
circularly polarized electromagnetic wave propagating in a plasma is investigated. A model based 
on the relativistic Boltzmann-Vlasov equation with a source term corresponding to the Schwinger 
formula for the pair creation rate is used. The damping of the wave, the nonlinear up-shift of its 
frequency due to the plasma density increase and the effect of the damping on the wave polarization 
and on the background plasma acceleration are investigated as a function of the wave amplitude. 
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The production of electron-positron pairs under the action of electromagnetic field (Schwinger effect) attracts great 
attention because it is a nonlinear effect that lies beyond the limits of perturbation theory. The study of this effect 
can shed light on the nonlinear properties of the quantum electrodynamics (QED) vacuum. 

This effect was first predicted for the case of a constant electric field more than 60 years ago (see also It is 

well known that a plane electromagnetic wave cannot produce electron-positron pairs because both its electromagnetic 
invariants, (E 2 — B 2 )/2 and E • B, are equal to zero. For this reason this effect was first considered in the case of a 
constant electric field, in which case the first invariant (E 2 — B 2 )/2 does not vanish. Later, this analysis was extended 
to the case of a spatially homogeneous time- varying electric field 0, IE IE 0> @ j but these results were long believed to 
~l ' be of academic interest only, because the power of the laser systems available at that time was far below the limit for 
pair production to become experimentally observable |o. 13. Il0| . However the recent development of laser technology 
has resulted in the increase of the power of optical and infrared lasers by many orders of magnitude . Presently, 
\ lasers systems are available that can deliver pulses with intensities of the order of 10 22 W/cm 2 in the focal spot. 
II ■ Such intensities are still much smaller than the characteristic intensity for pair production Isch — 4.6 x 10 29 W/cm 2 , 
which corresponds, for a laser pulse with wavelength « 1/im, to an electric field equal to the critical Schwinger field 
Esch = 1-32 x 10 16 V/cm. Nevertheless, there are projects that aim to reach intensities as high as 10 26 — 10 28 
W/cm 2 already in the coming decade. In addition, several methods for reaching the critical intensity with presently 
J> . available systems have been proposed recently. One of these schemes was demonstrated in the experiments at SLAC 
where 10 18 W/cm 2 laser photons, back-scattered by a 46.6 GeV electron beam, interacted with the laser pulse and 
several electron-positron pairs were detected [l^. Another scheme for reaching critical intensities was suggested 
in Ref . , where the interaction of the laser pulse with electron density modulations in a plasma, produced by a 
counterpropagating breaking wake plasma wave, results in the frequency up-shift and pulse focusing. In this scheme 
intensities of the order of the critical density can be obtained using 10 18 W/cm 2 laser pulses. Hence, a more detailed 
study of the Schwinger effect in time-varying electromagnetic fields and of all the processes that accompany it has 
become an urgent physical problem from an experimental point of view also. 

The process of electron-positron pair production by electro mag netic fields that are solutions of the Maxwell equations 
in a plasma and in vacuum was studied recently in Refs. |l4l Il5j | . In Ref. |15| it was shown that already for intensities 
of the laser pulse smaller than the critical intensity, the energy loss due to pair production is of the same order of the 
energy storied in the pulse. Therefore it is no longer possible to consider the electromagnetic field in the pulse as an 
external field and the energy loss by the electromagnetic field due to pair production and particle acceleration must 
be taken into account. 

The problem of the backreaction of the produced particles on the background field was discussed extensively in a 
number of papers on the particle formation process in high energy hadronic interactions 0, 0, 0, 0] as well as 
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under the action of electric fields |2]| . In the former case this process can be viewed as the quantum tunneling of 
quark-antiquark and gluon pairs in the presence of the background color-electric field of quantum chromodynamics 
(QCD). Such a color field is formed between two receding nuclei which are color charged by the exchange of soft 
gluons at the time of collision. This leads to the formation of a very strong color-electric field and, hence, to more 
copious pair production. It was understood that, in solving a dynamical problem with a strong initial electric field, 
the effect of the produced particles on the electric field (the back reaction) should be taken into consideration. The 
quark-gluon plasma emerging through tunneling in nucleus-nucleus collisions will change the color-electric field due 
to the appearance of conduction and polarization currents. The first current is due to the particle motion in the field 
while the latter arises from the process of pair creation. A kinetic equation coupled to Maxwell equations was used to 
solve this problem. The pair production process was considered as if occurring in QED, and the color-electric field was 
assumed to be Abelian, spatially homogeneous and time-dependent. However a spatially homogeneous time dep endent 
electric field is not a solution of Maxwell equations in vacuum. We also note that in Refs. [l^ll7LIT8LIT^.l20ll2l| special 
attention was paid to the properties of the emerging plasma, while the properties of the background field were not 
studied in detail. 

In the present paper, we consider the process of electron-positron pair production in a cold collisionless plasma, 
under the action of an electromagnetic field which is an actual solution of the Maxwell equations, as well as the 
backreaction of the produced pairs on the background field. In doing so we use the Boltzmann-Vlasov equation, with 
a source term obtained from the pair production rate 0, 0, ^| . In order to elucidate the role of the magnetic field 
component on the electron-positron pair production, we consider a planar, circularly polarized, electromagnetic wave 
propagating in an underdense collisionless plasma (for the sake of simplicity we consider an electron-positron plasma) . 
In the case of a plane wave in a plasma the first field invariant (E 2 — B 2 )/2 is not equal to zero due to the different 
dispersion equation with respect to that in vacuum. Therefore, in a plasma, electron-positron pairs can be produced 
by a plane electromagnetic wave, as was shown in Ref.|l4|. In this case a Lorentz transformation to the reference 
frame moving with the group velocity v g of the wave transforms the electromagnetic field into a purely electric field, 
that rotates with constant frequency, and with no associated magnetic field. Although this transformation reduces 
the problem under consideration to the situation where the pairs are produced by a time-varying electric field, the 
effects of the wave magnetic field are incorporated rigorously into our model. We notice that a similar approach was 
used earlier in Rcf . 22] . However in this latter paper a linearly polarized electromagnetic wave was considered and an 
approximation was adopted that is only valid in the limit of a small amplitude electromagnetic field. 

In the present paper we consider the effect on the background wave in the plasma caused by the pairs produced by 
the wave through their polarization and conduction currents. In particular, by considering the interaction between the 
wave and the plasma as an initial value problem in the moving frame, we study the evolution of the wave electromagnetic 
field. Due to the nonlinear properties of the equations governing the field evolution in time, we find a strong nonlinear 
dependence of wave field properties on the wave initial amplitude. We find a nonlinear up-shift of the wave frequency, 
a change of its polarization state and damping of its amplitude. In order to exemplify these effects, we consider two 
limiting cases of the electric field evolution. In the first limit the amplitude of the initial electric field is assumed to be 
so large that most of the electron-positron pairs are produced instantaneously. This leads to an instantaneous change 
of the electric field amplitude and frequency. In the second limit we consider a regime where the pair production rate 
is relatively small, so that all parameters of the wave change slowly. In this case we see the wave amplitude damp 
with time. 

This paper is organized as follows. In Sec^ we review some well known properties of a strong electromagnetic 
wave in a plasma. In Scc lIIII we study the equations governing the process of pair production and the evolution of the 
electromagnetic wave. In Secs lIVIVl we consider two limiting cases: fast changing field and slow changing field. The 
main results and conclusions are presented in Sec lVll 

II. RELATIVISTICALLY STRONG ELECTROMAGNETIC WAVE IN A PLASMA 

First we recover the properties of a relativistically strong electromagnetic wave propagating in an underdense 
collisionless electron-positron plasma that are needed in order to determine the form of the time varying electric 
field. This discussion is based on the results obtained by Akhiezer and Polovin in Ref. [23j. We consider a circularly 
polarized electromagnetic wave, propagating in an electron-positron plasma. In the following we use the c = 1 and 
h = 1 convention. 

An electromagnetic wave propagating in a plasma has a group velocity smaller than the speed of light in vacuum. 
Thus it is possible to make a transformation to the reference frame moving with the wave group velocity v g . In this 
frame the magnetic field of the wave vanishes and its time- varying electric field is spatially homogeneous and is thus 
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governed by the equation 



rfE 
~dt 



= -J?rj = -4,7 > J??- J v Q / Q (p,t)d 3 p. 
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where v = p/(m + £>)', p = |p| = (pi + P y + P z ) , fa(p,t) is the positron (electron) distribution function, 
normalized such that J" f a (p,t)d 3 p/(2ir) 3 = n a gives the number n a of electrons or positrons per unit volume, and 
e a is their electric charge with a = + for the positrons and a = — for the electrons. We consider the case of an 
electrically quasineutral plasma, where the density of positrons is equal to the density of electrons, n + = n_ = no, 
with no the density in the moving frame. If we assume that the plasma is cold, i. e., that the particle distribution 
function of the species a can be written as f a (p,t) = n a (2Tr) 3 S(p — p a (t)), the system of equations for the electric 
field evolution in the moving frame reduces to 
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We assume symmetric initial conditions, so that the positron and the electron momenta p±,± perpendicular to the 
direction of propagation of the laser pulse have opposite signs and equal absolute value, p+,j_ = — p = p_i_, while 
they have equal parallel momentum p|| = P||,o- Here the non zero value p\\^ = —mv g "f g of the parallel momentum is 
due to the Lorentz transformation from the laboratory to the moving frame and j g = (1 — v^/c 2 )^ 1 / 2 . Then, we find 
for the particle momentum p the equation 



fpi 

dt 2 



87re nop_L 
(m 2 +p 2 ) 1 /2 



The solution of this equation is 



p^ = — P(e x cos fit + e y sin fit), and A = Aq (e x cos fit + e y sin fit) , 



where A is the vector-potential, Aq = P/e and 



CI = 



87re 2 no 
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1/2 
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is the Langmuir frequency which enters the dispersion equation of an electromagnetic wave propagating in a plasma. 
The wave is propagating along z-axis. In this moving frame the wave electric field is given by 



E = flAo (e x sin fit — e y cos fit) . 
In the laboratory frame the wave dispersion equation is given by 

uj 2 = k 2 + n 2 . 



(7) 



(8) 



where f2 can be re-expressed in terms of the electron and positron densities n — n^/^ g and of the particle energy 
(r^ + P 2 ) 1 / 2 in the laboratory frame as ft = [(87re 2 n)/(m 2 + P 2 ) 1 / 2 ] 1 / 2 . The phase and group velocities of a nonlinear 
electromagnetic wave in a plasma depend on the plasma parameters and on the wave amplitude. From Eq.JSJ we find 
that the phase velocity v p h — u>/k and the group velocity v g = duj/dk, are related according to equation v p hv g = 1. 
In the laboratory frame the electric and the magnetic fields are given by 

E = —dtA = uAo[e x sm(u)t' — kx) — e y cos(ojt' — kx)}, (9) 



B = V x A = kA [e x cos(wi' — kx) — e y sm(u)t' — kx)], 
respectively, with t' the time in the laboratory frame. 



(10) 
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We see that in a plasma the first invariant of the electromagnetic field T = (E 2 — B 2 )/2 is not equal to zero and is 
given by 

rt 2 a2 i fn\ 2 2 , . 

It vanishes when the plasma density tends to zero, i. e., in vacuum. In the following we shall use the notation 
E = QA EE (n/uj)E . 



III. KINETIC DESCRIPTION OF THE ELECTRON-POSITRON PLASMA 



We consider the propagation of a circularly polarized electromagnetic wave in an underdense collisionless plasma 
in the reference frame moving with the wave group velocity v g . The relativistic kinetic equation 

^ + e a V^= qa (E,p), (12) 

describes the dependence on time and momentum of the distribution function f a (p,t) in this moving frame where 
a spatially homogeneous electric field E is present. The source term in Ea. ljT2*)l is proportional to the quasiclassical 
probability 



exp 



7r(TO 2 +P 2 L ) 



|eE(t)| 



(13) 



of tunneling through the gap between the lower and the upper continuum of electron energy spectrum in the presence 
of the electric field. We note that the form of Eq. l|13fl corresponds to the case of constant field 2]. However, the 
characteristic time of pair production c/l c , where l c = h/mc is electron Compton wavelength, is negligible with respect 
to the wave period. This estimate gives a lower bound on the pair production time while the estimate that follows 
from the quasiclassical approximation gives for the time of sub-barrier motion i tun = X/aio |?l.ll4|. where a = eA/mc 
is the dimensionless amplitude of vector-potential. However, for a 3> 1 (the case we are considering), even this 
estimate yields a pair production time much shorter than the wave period. Thus, it is possible to use Ea. (|13fl for 
the time-varying electric field with time playing the role of a parameter. In addition, following the reasoning of Refs. 
|l6l Il7l ITsj . we assume that the pairs are produced at Rest, i. e. the momentum distribution of the source term is 
taken to be proportional to the Dirac delta function 



q a {E,p) = 2e 2 E 2 (i)exp 



|eE(t)| 



S(p). (14) 



This assumption is reinforced by the fact that the momentum distribution in Ea. ^l3|) has a width p± ~ (| eE(t) | ) 1 / 2 = 
?7T,(|e(t) l) 1 ^ 2 <C m which is negligible with respect to the momentum that electrons (positrons) acquire in the elec- 
tric field. Here e = eE/m 2 = E/Es c h <gc 1 is the normalized electric field and E^ = m 2 /e is the critical 
Schwinger field. The source term has been normalized in such way that Jq a (E,p)d 3 p/(2tt) 3 gives the number 
[{eE) 2 /4-7T 3 ] exp[— irm 2 /eE\ of positrons (electrons) produced according to Schwinger's formula. 

We solve Ea. (|12l) by integrating it along the particle characteristics. The equations for the characteristics for each 
species a are 

dp± df a 

pii=pii,0) ~~dT = 6a ' ~dJ = qc " 

t 

Introducing the function A(t) = — j Eds, we obtain 

o 

t 

p_L = e a J Eds + pxo i-c, p± + e Q A(t) = p± . (16) 
o 

As a result the distribution function is given by the following expression 

t 

fa = f a ,o\p\\,Px + e a A(t)] + J q a {p x + e a [A(t)-A(t%t'}dt', (17) 
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where f a .o(p\\\,P±) is the distribution function of the initial plasma positrons (electrons) before the passage of the 
electromagnetic wave. Let us assume that at the initial time t = the plasma is cold so that 

f a , a =n (2ir) 3 5( P± )6( Pll -p l0 ), (18) 

where we recall that p\\ is the component of the particle momentum parallel to the electromagnetic wave propagation 
and pii o is its initial value which arises from the Lorentz transformation from the laboratory to the moving frame. 

The modification of the kinetic equation given by the source term in Ea. (|14|l must also be accompanied by a 
change of the source term in Maxwell equations. The electron-positron pair production by a spatially homogeneous 
time-varying electric field leads to the appearance of a time-dependent electric dipole which generates a polarization 
current. Thus the current density in Eq.Q acquires an additional term with respect to the situation when no pair 
production is present. Then Eq.|3J) reads as [l8| 



dt 

where the conduction current is 



dE 

— = -47Tj tot = -4tt (j cond + J p oi) , (19) 



d 3 p 

(m 2 +p 2 ) 1/2 (2tt) 3 ' 



E lf*M —* ■ _,m/, ( 2 °) 



and the polarization current is |16| (see the Appendix for details) 

a=+, 

Using the distribution function (|17|l . we obtain the following expressions for the current densities 

t 



= ? / qa{pMm2 +p2)1/2 w?- (21) 



A(t) „ , f Aft)-A(t') |eE(t') 
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'M2 



|eE(t')| 



dt', (22) 



2tt 



: E(t) exp 



'E(t)|. 



(23) 



where A = |A(t)|. In performing the momentum space integration we have used the fact that the pairs are produced 
with zero momentum. Inserting these expressions for the current densities into the r.h.s. of Ea. (|19fl and using the 
dimensionless vector-potential a = eA/m and the normalized electric field e = eE/m 2 , we obtain the equation for 
the electric field evolution in the presence of pair production 



( da(t) _ 
dt ~ 

deft) 

71 — ; — 

dt 



-me(t), 



a(t) 



p [1 + p 



em 
2^ 



eft) exp 
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a(t)-a(f) 



|e(t')| 2 



[l + |a(t) -a(t')| 2 ] 1/2 8tt 3 



■ exp 



|e(t')| 



|e(t)| 



dt' 



(24) 



1/2 

where w p = (87re 2 no/m) is the non-relativistic Langmuir frequency piiQ = p\\o/tn and k = 8ire 2 m 4 , where the factor 
m 4 stands for the inverse of the invariant Compton 4-volume m 4 = c/Z 4 w 0.14 x 10 53 cm -3 s _1 . Similar equations 
were obtained in Ref.|l8|. where however there was no initial distribution function and a spatially homogeneous 
electric field in vacuum was used, which is not a solution of Maxwell's equations. 

The nonlinear integro-differential equation l|24Jl cannot be solved analytically. Numerical solutions of this equation 
are presented in Fig.l for different initial amplitudes. We can see that the process of electron-positron pair production 
leads to the damping of the wave in the plasma and to the nonlinear up-shift of its frequency. The damping is due 
to the fact that each pair creation takes a portion of the field energy equal to 2mc 2 as well as the amount needed 
for the particle acceleration. The up-shift of the field frequency is due to the increase of the plasma density, and 
thus of the Langmuir frequency as new pairs are created. This frequency up-shift is seen in Fig.l, and bears a 
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FIG. 1: Time evolution in the moving frame of the x and the ^/-components of the dimensionless vector potential for different 
initial amplitudes: a — 1.4 x 10 5 (a), a — 1.5 x 10 5 (b), a = 1.9 x 10 5 (c) with initial plasma density no = 10 19 cm -3 in the 
moving frame; v g w 1, -y g = 10. The upper row shows the x-component of the vector-potential, the lower the y-component. On 
the X-axis time is measured in seconds; a = 1 corresponds, for a 1 /im wavelength pulse, to an intensity of 10 18 W/cm 2 and 
a — 4.6 x 10 5 to the Schwinger intensity. 



strong resemblance to the blue-shift of an electromagnetic wave that propagates in a medium that becomes ionized, 
as investigated theoretically in Ref.[2(| and experimentally in Ref . [2?!] . 

Since the pair production rate depends on the field amplitude exponentially, an unbalanced damping of the field 
components can occur and lead to a change of the field polarization. We should also note that the decrease of the 
amplitude of the vector-potential is accompanied by a frequency up-shift, so that the decrease of the amplitude of 
the electric field is not as fast as that of the vector-potential. These properties of the electric field are shown in 
Fig. 2, where the projections of the polarization vector are presented for the same set of initial parameters as in Fig.l. 
In Fig. 2a we see the damping of the x-component of the electric field and the transition from circular to elliptic 
polarization with the major axis of the ellipse directed along the y-axis. In addition, in Fig. 2b we see a rotation of 
the principal axes of the ellipse. The rotation of the principal axes of the polarization ellipse was discussed in Ref.[24| 
in the case of the free propagation (without source terms) of an elliptically polarized non liner pulse in a plasma. 
The situation shown in Fig. 2c is different from the two previous ones. In this latter case the pair production rate at 
the beginning of the field evolution is so large (see Fig. 3c) that the first wave oscillation cycle cannot be completed, 
leading to oscillations of the x-component of the wave vector potential around a non-zero mean value determined by 
the balance between the time averaged parts of the first two terms on the r.h.s. of the second of Eqs. (|24ll . This 
shift of the center of the oscillations of the x-component of the vector potential leads to a reduction of the oscillation 
frequency of this wave component so that, in this case, the x and the y-components of the wave oscillate at different 
frequencies. 

The increase with time of the pair plasma density is shown in Fig. 3 starting from an initial density before the 
start of the pulse propagation. We see that the particle density increases in steps as a function of time. This time 
dependence is due to the oscillations of the amplitude of the electric field that occurs in an elliptically polarized wave 
(see Fig. 2). The electron-positron pairs are mainly produced near the maxima of the electric field amplitude, while 
there is practically no pair production in between. From these plots we can also see that the frequency is up-shifted 
since the width of the steps decreases with time. For the pulse amplitude of frame c, the pair production occurs 
almost instantaneously at the beginning of the pulse evolution. 

The difference between the above three cases is clearly illustrated by the different shapes of the particle distribution 
functions in the p x -p y plane (Note that the electron and the positron distributions are one the mirror image of the 
other). In cases a) and b) electrons and positrons are mostly created at the maxima of the electric field |E| (and thus 
of the vector potential |A|). Since at birth = 0, in the case of a circularly polarized electric field this should lead 
to a ring type distribution. However, since the wave polarization becomes elliptical because of the backreaction due to 
the pair creation, the distribution function of each population consists, in the canonical momentum p± + e a A plane, 
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mainly of two blobs at ±e a A max . In the p x ~Py plane, these blobs move according to the time evolution of the vector 
potential A . On the contrary the position of the initial distribution function (denoted by a dark dot in the figure) 
corresponds to + e a A = 0. In case c) the pairs are created mostly at the start at p^ + e a A = e a A(t = 0). Since 
the time evolution of A(i) is ergodic, as shown in Fig. 2, their distribution tends to be randomized in the p x -p y plane. 



Py 
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FIG. 4: The electron distribution functions versus p x and p y in the moving reference frame for the same set of initial conditions 
as in Fig. 1 at time 2 x 10 -10 s. Particle momenta are normalized on the dimensionless vector-potential amplitude a multiplied 
times 10 5 . Black circles correspond to the initial plasma particle distribution at time 2 x 10 -10 s. 



The particle distribution function is shown in Fig. 5 versus the parallel momentum in the laboratory frame 
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a) b) c) 

FIG. 5: The electron distribution functions versus py in the laboratory reference frame for the same set of initial conditions as 
in Fig. 1 at time 2 x 10~ 10 s. Momentum is measured in the same units as in Fig. 4. The vertical lines correspond to the 
distribution of the initial plasma particles. 



at time t — 2 x 1CP 10 s. Note that in case c) the strong damping of the wave due to the pair creation and the 
resulting non adiabatic interaction has lead to a strong acceleration of the particles in the initial plasma. Such 
large values of the longitudinal momentum of electrons (positrons) in the laboratory frame are due to the transverse 
acceleration of electrons (positrons)in the moving frame. Performing the Lorentz transformation back to the laboratory 
frame we obtain for the longitudinal momentum in the laboratory frame of the initial electrons and positrons py = 
lg\P\\o + v g (l +P\ + |a 2 |) 1/2 ] « lgv g \a\ where we used |a| > |p|| |. 

In summary, the production of electron-positron pairs by the electromagnetic wave propagation in the plasma leads 
to the up-shifting of the wave frequency and to the damping of the wave amplitude and changes the polarization 
state of the wave. In order to illustrate these effects analytically, we consider two limiting cases of the electric field 
evolution, which can be clearly distinguished on the basis of the results shown in Fig. 1. First, we will consider the 
case of a fast changing field when the initial electric field amplitude is so large that the most of the electron-positron 
pairs is produced instantly. This case can be illustrated by the results shown in Fig. lc. Then, we will consider the 
case of a slowly changing field, when the pair production rate is relatively small and all the parameters of the wave 
change slowly as illustrated in Fig. la. 



IV. FAST CHANGING ELECTROMAGNETIC FIELD 



We assume that the amplitude of the initial electric field is so large that the electron-positron pair production and 
the consequent change of the properties of the wave electric field occur almost instantaneously. In this case we may 
approximate the time dependence of the source term as 



Qa(E,p,t) = Td(t)q a (E ,p), 



(25) 



where Eq is the initial amplitude of the electric field and r is a characteristic time defined so that the total number 
N of pairs produced is kept constant and given by r(|eo| 2 /47T 3 ) exp [— 7r/|eo|] per Compton 3-volume l^. Then, it is 
possible to carry out the integration over t' in the r.h.s. of Ea. (|24|l . As a result we obtain 



fa. 



a(t) 



■Pfo- 



>(t)]V2 



|e r 
8tt 3 



exp 



7T 

R 



a(i) - a 



[l + |a(i)-a | 2 ]V2 



(26) 



We note that ao in Ea. i|2t)[l and the initial value a(0) of the amplitude of the vector potential in Ea. H26|l are 
connected in an indirect way since the initial conditions for Ea. (|26|) are determined by the value of the electric field 
after the main part of the pair production has already taken place when the amplitude of the field has been reduced 
so that it is no longer capable of producing a significant amount of pairs. Therefore the values ao, eo that enter the 
r.h.s. of Ea. (|26[) should be considered as referring to the vector potential and to the electric field at the instant when 
the pairs are created. From these values we can deduce the value of a(0) using energy conservation and taking into 
account the polarization current. 

In order to understand the basic properties of the electric field evolution described by Ea. (|26() . it is convenient to 
linearize it. To do so, we assume that a <C 1 and obtain 



'(f) + (wi 2 + «JV)a(t) - uNao = 0, 



(27) 
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where uj' 2 = + p 2 ^) 1 / 2 . We see that Ea. (|27|l describes oscillations with frequency {uj' 2 + kN) 1 / 2 . Taking 

ap = (ao z , 0) we obtain 



a x (t) 
a y (t) 



CLQxKN 



CLQxKN 



to' 2 + kN uj' 2 + kN 



(lj> 2 + kN) 1 'H 



a\ cos 



—a\ sin 



(28) 
(29) 



where the initial conditions are a(0) = (oi,0), a'(0) = (0, {uj' 2 + kN) 1 / 2 ^). Thus, Eqs. QgEg) describe an elliptically 
polarized field oscillating with an upshifted frequency, as consistent with the numerical solution of Ea. i|24[l . 



V. SLOWLY CHANGING ELECTROMAGNETIC FIELD 



We assume that the pair production rate is small, so that the parameters of the electromagnetic wave change 
slowly. Then, we can neglect a(t) with respect to a(t'), since, as the wave amplitude decreases with time, a{t) 
becomes increasingly small compared to a(t'). Then, after differentiation with respect to time, Ea. l24fl reads 



d 3 a(t) 
dt 3 



a(t) 



a(t) 



p dt{l+p 2 +a 2 (i))V2 



»(t)]V5 



\e(t)\ 2 exp 



|e(t)| 



(30) 



If we linearize Eq. I|30[l and set the pair production rate equal to constant, we obtain 



a"'(t) 



n 2 V/'- 
P\\a) 



r a'(t) - «Wa(t) = 0, 



(31) 



where W = \e\ 2 exp (— 7r/|e|) and e is the electric field amplitude averaged over the field evolution. The equations for 
the two components of the vector-potential are decoupled and thus the polarization state is preserved. The solution 
of Eq. P|) is 



,{t) = C\ x y exp(— iwit) + C 2x , y exp(-icj 2 t) + C 3XiV exp(-zai 3 t), 



(32) 



where the frequencies Wj (i = 1, 2, 3) are the roots of the third order polynomial equation y 3 —yuj 2 /(l+p 2 Q ) 1 ^ 2 +iKW = 
and the constants Ci X>y are determined by the initial conditions. It is obvious from the form of the polynomial equation 
that one of its roots is imaginary and positive, while the other two are complex. The two terms in the solution with 
complex frequencies describe a damped wave, while the third term corresponds to a spurious, exponentially growing, 
term that can be excluded by an appropriate choice of the initial conditions. 



VI. CONCLUSION AND DISCUSSIONS 



In the present paper we considered the problem of the backreaction of the produced electron-positron pairs on the 
electromagnetic wave. We showed that there is a loss of wave energy due to the pair production and acceleration of 
these pairs in the electromagnetic field of the wave. 

We studied the propagation of a relativistically strong electromagnetic wave in an underdense electron-positron 
plasma. As is well known, a plane wave does not produce electron-positron pairs in vacuum. However the situation 
changes in a plasma. Due to the fact that in plasma the first field invariant T does not vanish, a plane wave can 
produce pairs that backreact on the wave. In order to describe the behaviour of the electrons and positrons, we 
adopted a kinetic plasma description and used the relativistic Boltzmann-Vlasov equation. Solving this equation we 
obtained the distribution functions of electrons and positrons which we used to derive the current density that enters 
the r.h.s. of Maxwell equation for the electric field. We note that all the calculations were carried out in the reference 
frame moving with the group velocity of the wave. In this frame there is no magnetic field as the wave has only an 
electric field. 

The solutions of the Maxwell equation for the evolution of the electric field in the plasma for different initial 
amplitudes are shown in Fig. 1. We see that, when the process of electron-positron pair production is taken into 
account, the evolution of the vector-potential in the plasma leads to the damping of the wave, accompanied by a 
nonlinear shift of its frequency. The damping of the vector-potential is due to the pair production which takes away 
a portion of the wave energy. This process is followed by the increase of plasma density which leads to the up-shifting 
of the wave frequency. 
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The damping of the vector-potential amplitude, together with its frequency increase, gives rise to an interesting 
phenomenon - the decrease of the wave amplitude followed by the increase of its frequency. Pair production decreases 
the wave amplitude while the growth of plasma density all over the space causes the increase of the wave frequency 
which is proportional to square root of the density. This effect resembles the case when the electromagnetic wave 
propagates in an ionizing medium |26|. 

Since the pair production rate depends on the field amplitude exponentially, an unbalanced damping of the field 
components occurs and leads to the change of the wave polarization (see Fig. 2). One can clearly see the change of 
the electric field polarization from circular to elliptic and the decrease of one of the two electric field components. 

Finally, we considered two limiting cases in order to identify the properties and the exact causes of this nonlinear 
behaviour of the wave. In the first case the initial amplitude of the electric field is so large that the pair production 
and the change of the wave properties occur instantaneously. In the second case the pair production rate is small so 
that the properties of the wave change slowly. In first case we obtain an electric field oscillating with a new frequency 
and amplitude. The frequency increase is due to the fact that the density of plasma increases because of the pair 
production. In the second case we obtain a damped wave with an amplitude that decreases slowly with time. We see 
that the process of electron-positron pair production leads to the damping of the wave. There is a nonlinear shift of 
its frequency due to the decrease of its amplitude: as a result the frequency is up-shifted. 
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Appendix: Properties of the polarization current 



The form of the polarization current can be derived from the energy-conservation law. In order to calculate the 
expression for the polarization current we calculate the second moment of the kinetic equation (|12fl . 



2 i 2\l/2 

m + p ) ' — 



a=± 

which we write as 



where 



9 fa tp^/q i ,.; 



dK 



a— ± 



2 i 2\l/2 

m + p ) 1 



q a (E,p)d 3 p, 



. dn 

-57 — JconrfE = £ — 2m—, 



dt 



at 



(33) 



(34) 



K =H /[(m 2 +p 2 ) 1/2 -ml/ a (p,t)d 3 p 

ot=± J 



and 



Jcond ^ 



/ (m2+ p p2)1 / 2 E e «/«(p»*)^ 



is the conduction current. The second moment of the r.h.s. of kinetic equation r.h.s. gives two terms. One of them, 
—2mdn/dt, is related to the rest energy increase due to the pair production, and the another is equal to 



£=E [(™ 2 +p 2 ) 1/2 q Q (E,p)d 3 p. 
a=± J 



We can thus rewrite Eq. (|34|l in the following form 

d 



dt 



(K + 2nm) = j conrf E + £, 



(35) 
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where the expression inside the brackets correspond to the full energy of the plasma particles. The energy balance 
equation for the electromagnetic field in the moving reference, frame where only the electric field is present, can be 
obtained by multiplying Eq. (fP3|) by the vector E 

d /E 2 \ 

dt \ 8vr ) = ~' icondE ~ j po ' E ' ^ 36 ^ 
Here we took into account that no external current is present in our case. Adding Eqs. (|35|) and l|3t)|) we obtain 

^ (k + 2n + = -J P o/E + S. (37) 

The expression inside the brackets in Eq. (|37|l is the energy of the system of the electromagnetic field and the plasma 
electrons and positrons. Energy conservation requires the r.h.s. of Ea. l|37|) to be zero. Thus the polarization current 
should be of the form 



E I 

hoi = t^E / (™ 2 + p 2 ) 1/2 fe(p, t)d 3 p. 



(38) 
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